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Abstract. We explore a simple model for the representation of the observed distributions of the 
motions, and the characteristic ages of the local population of pulsars. The principal difference from 
earlier models is the introduction of a unique value, S, for the kick velocity with which pulsars are 
born. We consider separately the proper motion components in galactic longitude and latitude, and 
find that the distributions of the velocity components parallel and perpendicular to the galactic plane 
are represented satisfactorily hy S = 200 km/sec, and leave no room for a significant fraction of much 
higher velocities. The successful proposition of a unique value for the kick velocity may provide an 
interesting tool in attempts to understand the physical process leading to the expulsion of the neutron 
star. 

1 Introduction; the Model used 

The early proper motion measurements have indicated that pulsars, on the average, have spatial 
velocities of the order of a few hundred km/sec (Lyne, Anderson & Salter 1982). Subsequent mea- 
surements, over the past 15 years have provided evidence for speeds as high as a thousand km/sec 
(Bailes et al. 1990; Fomalont et al. 1992; Harrison et al. 1993). Many models have been suggested to 
explain the spatial velocities of pulsars. Harrison & Tadimaru (1975) proposed the "rocket" theory 
which essentially stated that the pulsar is accelerated along the magnetic axis due to the radiation re- 
action soon after it was born. However, subsequent observations have failed to confirm the prediction 
of the model, that there is an alignment between the directions of proper motion and the magnetic 
axis. Gott, Gunn & Ostriker (1970) proposed that pulsars derive their spatial velocities from their 
progenitor binary systems, when the binary disrupts due to the heavy mass loss during the supernova 
explosion (Blaauw 1961). Asymmetric supernova explosions were proposed by Shklowskii (1970). No 
fully satisfactory theory has been proposed yet. 

Clearly, knowledge of the shape of the birth-velocity distribution is of basic importance for any 
theory that attempts to clarify the mechanism that causes the birth-velocities. The present paper is 
concerned with this problem. Among the most recent hypotheses for the birth-speed distribution are 
those by Lyne & Lorimer (1993) and Hansen & Phinney (1997), shown in figure |l|. Hartman (1997) 
finds that a velocity distribution that has many more pulsars at low velocities than that of Lyne & 
Lorimer also describes the observations quite well. 

Whereas these distributions are characterised by maximum density in velocity-space near zero- 
velocity, we assume all pulsars to have been ejected from the parent body with the same birth speed 
- henceforth denoted by S' - isotropically distributed. Accordingly, in velocity space, velocity vector 
points define a spherical surface with radius 5'. We consider such an alternative worth exploring 
because it cannot be excluded at this stage of pulsar physics, that the birth (kick) velocities are 
spread around a (physically) preferred value rather different from zero. With regard to the space 
distribution of the progenitors, we assume these to be located in a layer of zero thickness at z = 0. 
This assumption is rather harmless because it is generally agreed that pulsars originate from the 
massive stars in the Galaxy which are confined to a layer the thickness of which is very small as 
compared to the scale-height of the z-distribution of the pulsars. Finally we assume an upper limit 
to the age, T, and we do not consider for this analysis pulsars with ages greater than T. Accordingly, 

^This paper is an extended analysis of the talk presented by AB at the Raman Research Institute on 20th February, 
1996. 




Figure 1: The birth speed distribution of pulsars as predicted by Lyne & Lorimer (1993) (soHd Une) 
and Hansen & Phinney (1997) (dash Une), and as explored in the present paper (gray bar). 

the quantity to be solved for will be only S. Our analysis also differs from previous ones in that we 
separately consider the velocity components parallel and perpendicular to the galactic plane. 

Naturally, the assumed uniqueness of the birth speed cannot but be an approximation of the 
true situation, for even if there is uniqueness in the kick velocities with which the pulsars leave the 
parent body, these are vectorially added to the space motion of the parent body, i.e., in the case of 
components of binary systems, the orbital velocity and the systemic velocity. This implies that the 
birth-speed distribution is not a delta function but exhibits a certain width. However, since the birth 
speed we arrive at, S = 200 km/sec, generally dominates over these vector additions, the simplified 
model remains a good first approximation. 

In section ^ we give analytical expressions for the expected distributions, and demonstrate why 
the pulsar velocities parallel to the plabe provides a direct way to estimate the velocity S. The 
observed pulsar population suffers from severe selection effects, which must be taken into account 
when comparing theoretical and observed distributions. We have done this through a detailed Monte 
Carlo simulations described in section |^. 

We compare our model with observations for the cylindrical volume of space defined by \z\ < 1 
kpc and projected distances dcosb smaller than 2 kpc only. Within this volume we assume steady 
state for velocity and space distribution. 

2 Definitions and analytical relations 

We denote by, 

5, the (unique) three dimensional speed at birth of the neutron star. In our analysis this will be 

the velocity with respect to the local standard of rest at the position of the Sun. 

T, the (adopted) upper limit to the pulsar age; only those pulsars with characteristic ages less than 

T are considered for this analysis. 

z, the distance from the galactic plane, 

K{z), the acceleration perpendicular to the plane due to the galactic gravitational field; K will be 
counted positive in the direction towards the plane, 
d, the distance from the Sun, 

w, the velocity component perpendicular to the galactic plane. 
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Table 1: The maximum height from the galactic plane (zmax) and the quarter oscillation time (T1/4) 
for various initial 2: -velocities (wq)- 



wq, the value of w at birth, 

V, the velocity component (parallel to the plane) corresponding to the proper motion component in 
galactic longitude, 

G{S), the number of neutron stars generated per unit time per imit surface at z = 0, 

Gq{wq) dwQ, the number among these in the interval wq, wq + dwQ. Since the velocities at birth 

are distributed isotropically, this is a flat distribution: 

GoK) dwo - G{S) dwo/S = C dw^, (1) 

Gz(w) dw, the number passing level z with velocities w, w + dw per unit surface per unit time. We 
assume symmetry with respect to the galactic plane, hence the model deals only with positive values 
of z and w. 

There is a subset of the pulsar population moving toward the plane after having reached maximum 
distance z. Table |l| gives this maximum distance, Zmax for various values of wq, and the quarter 
oscillation times T1/4, computed with the adopted force field described in section 3.1. Objects with 
Wq below about 52 km/sec remain during their full lifetime below z = 1 kpc and hence, within the 
volume of space we study. For S = 200 km/sec these represent about one quarter of the population. 
For T = 50 Myr, roughly one half of them will be observed also after they have reached maximum 
distance from the plane. Objects with wq in excess of 50 km/sec will orbit beyond z = 1 kpc and 
re-enter the lower domain on their way back. In these upper parts of their orbits they dwell relatively 
long: for wq = 60 km/sec they re-enter only at the age of about 50 Myr, for larger wq they dwell even 
longer beyond 1 kpc. Hence, it is only in the interval of wq between 52 and 60 km/sec that we expect 
to find objects younger than 50 Myr to re-appear in the domain below z — 1 kpc. Summarizing, for 
T = 50 Myr, a small fraction of the observed population, about 15 percent, will be on their way back 
to the plane. For lower values of T, this percentage is accordingly lower. The simulations take their 
presence into account. 

2.1 Velocities perpendicular to the plane 

The space density of neutron stars at z = with wq, + dw^ is 

Do{wo) dwo = (l/wo)G'o(wo) dwo = (C/wq) dwQ. (2) 

Since ^ 

w'^ = wl-2 I K{z) dz (3) 

for objects to reach level z, a minimum value of wq, denoted by Wmin is required: 

<n. = 2 r K{z)dz = 2P{z) 
Jo 

Where, P(z) = ( K{z) dz. (4) 

Jo 
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The range of values Wq for objects reaching level z and beyond is: 

Range in Wq = S'^ — w'^^^ (5) 
This is also the range of the corresponding values w'^ at level z, since Wmin becomes zero and S becomes 

where = - 2 [ K{z) dz = S"^ - 2P^ = - wf„i„ (6) 
Jo 

However, the range of w at level z, i.e., from zero to Ws, is broader than the corresponding range at 
birth in the ratio Wg/ {S — Wm). For instance, if we assume S=180 km/sec, for reaching level 420 pc, 
wq should range from 30 to 180 km/sec, i.e., 150 km/sec, whereas the resulting range at level z=420 
pc becomes w = to Wg — 177 km/sec. 

Since for given z, from equation |^, Aw'^ — Ao(wq), we have, wq Awq = w Aw. Hence, the 
broadening for given w is in the ratio wq/w. 

Accordingly, whereas G'o(wo) dw is flat, for a given level z, Gz(w) dw is proportional to w/wq, 
with w and wq related by equation ^. 

The number of objects per unit volume at level z with w in the interval w, w + dw will be, 

D,{w) dw = il/w)G,iw) dw. (7) 

Due to the broadening just mentioned, 

G,iw) = iw/wo)Goiwo) (8) 

Hence, Dz(w) = (C/wq), and we have the well known relation, Dz{w) — Dq[wq) for pairs of values w 
and Wq related by equation ||. 

The above relations do not take into account the disappearance from the sample of the pulsars 
older than age T. We denote by Wt the minimum velocity at z = required for reaching level z within 
the time-span T. wt is defined by: 

T = = f (9) 

io Jo ^Jwt-2P' ^ 

Thus, there is an interval of velocities at level z from zero to w^, in which no objects will occur. Again, 
we have, 

wl = wj - 2P 

Accordingly, at level z, the distribution of velocities w can be described by the above Dz{w) dw 
truncated at the low velocity end by eliminating the interval w — to w — Wz- 

Figures |^ demonstrate this effect. They have been obtained numerically for an adopted value 
S — 180 km/sec and T — 8 Myr. The five histograms correspond to five different z-slabs (between 
and 1 kpc), each of width 200 pc. In the highest layer we note the absence of pulsars in the velocity 
range w = to 90 km/sec. At lower z-levels the empty interval is accordingly smaller. At the high- 
velocity end the limiting velocity wi is fixed by wf — (S*^ — 2P). The fact that the cut-off is not sharp 
is due to the combination of all z-values within the slab. 



2.2 Velocities parallel to the plane 

We denote by, 

V, the component of the velocity parallel to the galactic plane corresponding to the proper motion in 
galactic longitude: 

V = d^i cos b 

Since the distribution of v is symmetric we consider only the absolute value of w. At level z, for a 
given T, v ranges from to vi = {S^ — w^Y^'^. The predicted distribution will be shown to be: 
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Figure 3: Distribution of velocities D{v) parallel to the plane, as given by equation ^ The three 
curves correspond to z = 0.1, 0.5 & 0.7 kpc, respectively, for S = 180 km/sec and T = 8 Myr. 



D{v)dv = C S'<^7r/2-arcsin 



Wt 



(10) 



where C is a constant defined below. 

We define by N{v, w) the number of objects at level z per unit of v and of w. Equation |l^ then 
follows from 

R dw 

D{v)dv ^ dv N{v,w) dw ^ dv N{0,w) 

Jwt Jwt 



Vi?2 - w2 



where, i?^ = {S^ - w^), wf = (S^ - v^), and N{0,w) = N{0,0) S/VS^-v^, so that. 



D{v)dv = dwS'iV(0,0) 



\/ ~ v"^ ~ w"^ 



S N{0,0) 



7r/2 — arcsin 



Wt 



dv. 



where C = iV(0, 0) is the number per unit v and unit w at t; = 0, ui = at z = 0. 

For low z-levels the distribution D{v) dv is nearly flat up to v close to vi. This is due to the fact 
that the second term in equation |lO| becomes significant with respect to the first, constant term, only 
at high V. For instance, for S — 200 km/sec and T = 20 Myr, a,t v = 0.9vi, the ratio between the 
second and the first term near z = is only 0.036. For high z-levels the distribution becomes less flat. 
Figure ^ demonstrates this effect. As we chose higher values of T, the nature of the curves for higher 
z becomes more similar to the one for the lowest z. The truncated nature of D{v), particularly at 
low z, provides an important tool in estimating S] the extreme observed value, vi is an approximate 
measure of 5" virtually independent of the assumed life time T. 



3 Simulations 



We did a Monte Carlo simulation to generate a large number of objects in the solar neighbourhood 
(with initial value of z—0), keeping the surface density constant. As mentioned in the previous 
sections, we assumed a simplistic model with all neutron stars being ejected out isotropically with 
a constant speed (S). This assumption corresponds to having initial velocities in the z-direction 
ranging from zero to S with uniform probability. We evolve them in an assumed local gravitational 
potential to get their velocity and spatial distributions, and in addition we studied the distributions of 
their characteristic ages. Each object was assigned an age chosen randomly with uniform probability 
between zero and T. Given the distribution of characteristic ages of known pulsars, we chose a value 
of T = 50 Myr as an optimal choice in order to include a sufficiently large sample, and at the same 
time describe the kinematics with our potential model. The presence of the objects which turn over 
after reaching their maximum z-amplitudc of oscillation within their age are accounted for. 



3.1 Local gravitational potential 

To study the kinematics of the generated objects we assumed the local gravitational potential function 
given by Kuijken & Gilmore (1989). The vertical acceleration due to this potential is given by 
(Bhattacharya et al.) 

= 1.04 X 10"^ 

where Qz is in units of kpc/Myr^ and z in kpc. We use this function to evolve every one of the 
generated objects from their initial position and velocity to their current value after a time t. 



1.26 z 



0.58 z 



(11) 



3.2 Further assumptions for simulation 

We assumed a Gaussian distribution of fields in logarithmic scale, with < log B{G) >— 12.2 with a 
root mean square spread around the mean of 0.3. The choice of initial periods will be discussed in 



section 4.4. With an assumption that the magnetic field does not decay significantly within the time 
scales of our interest (^50 Myr) the final rotational period of the pulsar can be easily calculated with 
the simple dipole formula 

2 ^ ^l^PP (12) 

where c is the velocity of light, / is the moment of inertia of the neutron star, and P is the time 
derivative of the rotation period P. 

The simulation was repeated for various values of S (ranging from 90 km/sec to 350 km/sec) 
and T (from 15 Myr to 50 Myr). At the end of the evolution, pulsars which have the quantity 
{B/P"^) < 2 X 10^^ were neglected from the simulation, since it is believed that the pulsar activity 
ceases to continue below this limit. 



3.3 Selection filter 

After evolving the generated objects in the galactic potential, their final position, velocity components, 
rotation period and characteristic age are noted down. After compensating for observational selection 
effects a subset of this population was selected to compare with the known sample. The procedure we 
have adopted for this task is the same as the one described by Deshpande et al. (1995). We selected 
only those pulsars which are, in principle, detected by any one of the eight major pulsar surveys 
considered, namely (1) Jodrell Bank survey, (2) U. Mass-Arecibo survey, (3) Second Molonglo survey, 
(4) U. Mass-NRAO survey, (5) Princeton-NRAO Phase I survey, (6) Princeton-NRAO Phase II 
survey, (7) Princeton-Arecibo survey, and (8) Jodrell Bank-1400 MHz survey. Since our analysis 
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deals with only the solar neighbourhood (projected distance onto the plane less than 2 kpc, and 
|2;| < 1 kpc) the effect due to interstellar scattering and dispersion are ignored. While comparing the 
results of the simulation with the known sample, only the subset of the simulated sample (detectable 
simulated sample) was used. 

For consistency, even the sample of known pulsars were passed through this filter to select a 
subsample, which was finally used for comparison with the simulated sample. 



4 Comparison with observations 
4.1 Velocities parallel to the plane 

As the first step, we consider the distribution of velocity components parallel to the plane, D{v) 



because as it was explained in section 2.2, these provide the estimate of S. From the compilation of 
Taylor et al. (1995) we use those proper motion measurements for which the quantity dxa < 94 km/sec 
(20 A.U./yr), where a is the formal error in the proper motion quoted by those authors. We do not 
apply corrections for solar motion and differential galactic rotation, since within the range of distances 
considered, these are negligible as compared to the true proper motions and their observational errors. 
Figure ||(a) gives the distribution of these velocities, and figure U(b) gives their distribution as a 
function of the height from the galactic plane. We note that with the exception of 5 objects, all 
are confined within ±210 km/sec. In comparison with figure ^ we may conclude that the value of 
the unique velocity must be roughly around about 200 km/sec. Clearly, values of S in excess of 250 
km/sec are incompatible with the proper motions. 



4.2 Density Distribution perpendicular to the plane 

Hartman & Verbunt (1994), from their study of evolution of neutron stars in various models of the 
local gravitational potential, demonstrate the effect of observational selection for understanding the 
distribution of pulsars perpendicular to the plane. Since observations tend to detect more luminous 
pulsars (which are younger and closer to the plane), irrespective of the nature of the true distribution 
the observable distribution is roughly the same. This has been seen in our simulations too. While 
comparing the simulated sample (after compensating for the selection effects) with the observed, a 
satisfactory Kolmogorov - Smirnov probability was achieved for many different combinations of S 
and T. For S = 180 km/sec and T = 50 Myr (figure |) the K-S probabihty was 59.6%, and for 
many other combinations of {S,T), like (250,50), (350,50) the K-S probability was roughly as high as 
for (180,50). Therefore, one can conclude that the density distribution perpendicular to the galactic 
plane alone is not a severe constraint on the kinematic properties of pulsars, and one needs to study 
other properties like the velocity components and the distribution of characteristic ages, etc. 



4.3 Velocities perpendicular to the plane 

We do not have measured w, since there is no way of measuring the radial component of the velocity 
vector. However, for pulsars with sufficiently low galactic latitude (say, |&| < 50°), we can approximate 
w — fib X d provided cos 6 is close to unity. Out of the 30 pulsars selected, 23 have cos 6 > 0.87, 5 
have cos b = 0.78 — 0.87, and 2 have cos b — 0.72 — 0.78. Figure ^(a) gives the distribution of [fn, x d). 
Figure ^(b) shows the distribution of /if, x d against the height from the galactic plane. Objects at 
negative z are plotted with reversed z and reversed /if, x d. The asymmetry in the distribution of 
fib X d in figure ^(a) shows the flow of objects away from the galactic plane. 

Once a considerably richer sample of accurate pulsar proper motions is available, it will be of great 
interest to make the detailed comparison between the observed distribution and the theoretical one 
of figure ^ and check its rapidly varying shape with increasing z. The predicted run of the mean 
velocity v with increasing z is in harmony with observations, as shown in figure @(b). 
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Figure 4: Velocities parallel to the plane corresponding to the proper motions in galactic longitude, (a) 
Distribution of these velocities, (6) Distribution of velocities as a function of height from the galactic 
plane. Only those pulsars whose proper motion errors {la) are less than 94 km/sec (20 A.U./yr) are 
considered for this figure. 
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Figure 5: The cumulative number distribution of pulsars as a function of height from the galactic 
plane. Dash line is for the known sample and the solid line is the best fit corresponding to S = 180 
km/sec and T = 50 Myr. 

4.4 Distribution of characteristic ages 

As described in section we have assumed that the magnetic fields at birth are distributed as a 
Gaussian in logi? around 12.2 with an r.m.s. of 0.3. Since many of the earlier works (Vivekanand 
& Narayan 1981; Srinivasan et al. 1984; Narayan 1987; Narayan & Ostriker 1990; Deshpande et 
al. 1995) suggest that pulsars are born with rotation periods as long as a few hundred milliseconds, 
for the simulation, all pulsars were assumed to be born with rotation period of 0.1 sec. With the 
assumption that the magnetic field does not decay with time, the period was evolved with the simple 
dipole formula |l^. Objects were evolved in the local gravitational potential for their corresponding 
ages (distributed between zero and T with uniform probability) . At the end of the simulation pulsars 
with B/ P'^ < 2 X 10^^ Gs^^ were neglected, since the pulsar activity is believed to cease below this 
limit. After applying the selection effect filter^ a subset of the simulated sample (detectable subset of 
the simulated sample) was selected for comparison with the known sample. 

Figure gives the distribution of known pulsars in the z — log Tch plane (where Teh is the charac- 
teristic age of the pulsar). The 'dash' line gives the observed mean value of log Tch in the five different 
z-slabs from zero to 1 kpc, and the solid line gives the corresponding mean value of the "detectable" 
simulated sample for the combination of S" = 180 km/sec and T = 50 Myr. The "dotted" line gives 
the possible lower limit of Tch for a given height from the plane, i.e., with initial velocities, wq, ranging 
from zero to 5, the object can reach a given height z, only if its age is greater than a certain value. 

As one can see from this figure, our model fits with the observed distribution of log Tch as a function 
of z, quite satisfactorily. The fact that the envelope in figure ^ (given by the dotted line) encloses all 
the pulsars means that almost all pulsars are born close to the plane and they migrate by the virtue 
of their velocities. As it turns out, to constrain the combination of S and T, the distribution of log Tch 
can be used effectively. 

Ideally we would have liked to compare the distribution of log Tch in each z-slab with the cor- 
responding observed distribution. However, small number of the observed sample prevents us from 
doing so. 
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Figure 6: (a) Distribution oidx jih of known pulsars, (6) distribution of known pulsars as a function 
of d X /if, and the height from the plane. Only those pulsars with |&| < 50° and velocity error less 
than 94 km/sec are included. The solid line gives the mean of the simulated sample for S = 180 
km/sec.,T = 50 Myr. 
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Figure 7: Distribution of known pulsars in z — log Teh plane. The dash line indicates the average 
log function of z, for the known sample. Solid line indicates the corresponding quantity for 

the simulated sample for S'=180 km/sec and T=50 Myr. The dotted lines give the predicted lower 
bound for function of z.The error bars indicate the uncertainty in the mean value. 

5 Conclusions 

We have explored the hypothesis that pulsars are born in the galactic plane with a unique kick 
velocity S, distributed isotropically. We therefore investigate the pulsars in the solar neighbourhood 
by studying the distributions of the velocity components parallel and perpendicular to the galactic 
plane, and the characteristic age distribution. We find that the known sample of pulsars in the solar 
neighbourhood is consistent with kick velocities of about 200 km/sec. 

We wish to emphasis that our analysis of the nearby pulsars leaves no room for the existence of an 
appreciable fraction of pulsars with velocities much higher than 200 km/sec. These should have shown 
up in the plots of figures ^ & ^. Evidence for very high velocities came from a few reliable proper 
motion measurements (Bailes et al. 1990; Harrison et al. 1993; Fomalont et al. 1992) and velocities 
inferred from the association of a few pulsars with supernova remnants (Caraveo 1993; Frail, Goss & 
Whiteoak 1994). The rehability of such associations has been questioned by Kaspi(1996). 

The analysis by Brandt & Podsiadlowski (1995) shows that the assumed kick speeds of the order 
of 450 km/sec are inconsistent with the eccentricity - orbital period distribution of binaries which 
contain neutron stars, and it is consistent with kick speeds of 200 km/sec. This supports our result. 

The velocity distribution derived by Lyne & Lorimer (1993) and Hansen & Phinney (1997), referred 
to in the introduction and figure 1, show a considerable fraction of pulsars with velocities around 200 
km/sec, the value of S we arrived at. They differ from our model in two respects: 

1. Both of them show maximum density in velocity space around velocity zero, whereas our model 
gives zero density at zero velocity. 

2. Both of them show a long tail towards the high-speed end. This is recognised in our plots of 
the velocity components v and w (figures 4 & 6), but this represents a small fraction only of the 
total distribution. 

The recent paper by Fomalont et al. (1997) has brought out improved proper motion measurements 
for about 15 pulsars. Among those pulsars with distance less than 2 kpc and galactic latitude less than 
50 degrees, on the basis of the earlier measurements we had ignored 1556-44 & 1749-28, since they had 
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errors greater than 20 A.U./yr. For 1919+21 & 2045-16 the improved proper motion measurements 

seem to be quite different from the earlier measurements. However, since these changes are \mhkely 
to change the results in this paper we have done our analysis on the basis of the earlier measurements. 

To summarise the main result, we show that if pulsars are postulated to derive a unique isotropic 
supernova kick speed at their birth around 200 km/sec, the distribution of the proper motions and 
the characteristic ages of pulsars in the solar neighbourhood can then be understood satisfactorily. 
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Figure 2: Distribution of velocities perpendicular to the plane of pulsars in different z-slabs after 
evolving in the potential for ages up to T. This plot corresponds to S = 180 km/sec and T = 8 Myr. 
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